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Ergodicity of the 

2D Navier-Stokes equations 

with degenerate stochastic forcing 

o 

O ' By Martin Hairer and Jonathan C. Mattingly 

^; Abstract 

< 

^D . The stochastic 2D Navier-Stokes equations on the torus driven by degen- 

^^ ' erate noise are studied. We characterize the smaUest closed invariant subspace 

for this model and show that the dynamics restricted to that subspace is er- 
Pm ' godic. In particular, our results yield a purely geometric characterization of 

•"^ . a class of noises for which the equation is ergodic in Lq(T^). Unlike previous 

■ ^ I works, this class is independent of the viscosity and the strength of the noise. 

2 ' The two main tools of our analysis are the asymptotic strong Feller property, 

introduced in this work, and an approximate integration by parts formula. The 
first, when combined with a weak type of irreducibility, is shown to ensure that 
cn I the dynamics is ergodic. The second is used to show that the first holds un- 

[»^ ' der a Hormander-type condition. This requires some interesting nonadapted 

OO ■ stochastic analysis. 

o 
^. 

O . 1. Introduction 

o 



In this article, we investigate the ergodic properties of the 2D Navier- 
Stokes equations. Recall that the Navier-Stokes equations describe the time 
evolution of an incompressible fluid and are given by 



-(— > 

s 

'^\ (1.1) dtu + {u-V)u = uAu-Vp + C , divu = 0, 

^\ ' where u{x, t) G R denotes the value of the velocity field at time t and position 

j^ ■ X, p{x,t) denotes the pressure, and S,{x,t) is an external force field acting on 

the fiuid. We will consider the case when x E T^, the two-dimensional torus. 
Our mathematical model for the driving force ^ is a Gaussian field which is 
white in time and colored in space. We are particularly interested in the case 
when only a few Fourier modes of ^ are nonzero, so that there is a well-defined 
"injection scale" L at which energy is pumped into the system. Remember 
that both the energy ||u|p = J \u{x)\'^ dx and the enstrophy ||V A n|p are 
invariant under the nonlinearity of the 2D Navier-Stokes equations (i.e. they 
are preserved by the flow of(l.l)ifi/ = and ^ = 0). 
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From a careful study of the nonlinearity (see e.g. [Ros02] for a survey 
and [FJMR02] for some mathematical results in this field), one expects the 
enstrophy to cascade down to smaller and smaller scales, until it reaches a 
"dissipative scale" r/ at which the viscous term i/Au dominates the nonlinearity 
{u-V)u in (1.1). This picture is complemented by that of an inverse cascade of 
the energy towards larger and larger scales, until it is dissipated by finite-size 
effects as it reaches scales of order one. The physically interesting range of 
parameters for (1.1), where one expects to see both cascades and where the 
behavior of the solutions is dominated by the nonlinearity, thus corresponds to 

(1.2) KL-i<r/~i. 

The main assumptions usually made in the physics literature when discussing 
the behavior of (1.1) in the turbulent regime are ergodicity and statistical 
translational invariance of the stationary state. We give a simple geometric 
characterization of a class of forcings for which (1.1) is ergodic, including a 
forcing that acts only on 4 degrees of freedom (2 Fourier modes). This charac- 
terization is independent of the viscosity and is shown to be sharp in a certain 
sense. In particular, it covers the range of parameters (1.2). Since we show that 
the invariant measure for (1.1) is unique, its translational invariance follows 
immediately from the translational invariance of the equations. 

From the mathematical point of view, the ergodic properties for infinite- 
dimensional systems are a field that has been intensely studied over the past 
two decades but is yet in its infancy compared to the corresponding theory 
for finite-dimensional systems. In particular, there is a gaping lack of results 
for truly hypoelliptic nonlinear systems, where the noise is transmitted to the 
relevant degrees of freedom only through the drift. The present article is an 
attempt to close this gap, at least for the particular case of the 2D Navier- 
Stokes equations. This particular case (and some closely related problems) 
has been an intense subject of study in recent years. However the results 
obtained so far require either a nondegenerate forcing on the "unstable" part 
of the equation [EMSOl], [KSOO], [BKLOl], [KSOl], [Mat02b], [BKL02], [Hai02], 
[MY02], or the strong Feller property to hold. The latter was obtained only 
when the forcing acts on an infinite number of modes [FM95], [Fer97], [EHOl], 
[MS05]. The former used a change of measure via Girsanov's theorem and the 
pathwise contractive properties of the dynamics to prove ergodicity. In all of 
these works, the noise was sufficiently nondegenerate to allow in a way for an 
adapted analysis (see Section 4.5 below for the meaning of "adapted" in this 
context). 

We give a fairly complete analysis of the conditions needed to ensure the 
ergodicity of the two dimensional Navier-Stokes equations. To do so, we em- 
ploy information on the structure of the nonlinearity from [EMOl] which was 
developed there to prove ergodicity of the finite dimensional Galerkin approx- 
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iniations under conditions on the forcing similar to this paper. However, our 
approach to the full PDE is necessarily different and informed by the pathwise 
contractive properties and high/low mode splitting explained in the stochas- 
tic setting in [Mat98], [Mat99] and the ideas of determining modes, inertial 
manifolds, and invariant subspaces in general from the deterministic PDE lit- 
erature (cf. [FP67], [CF88]). More directly, this paper builds on the use of 
the high/low splitting to prove ergodicity as first accomplished contempora- 
neously in [BKLOl], [EMSOl], [KSOO] in the "essentially elhptic" setting (see 
section 4.5). In particular, this paper is the culmination of a sequence of pa- 
pers by the authors and their collaborators [Mat98], [Mat99], [EHOl], [EMSOl], 
[Mat02b, Hai02], [Mat03] using these and related ideas to prove ergodicity. Yet, 
this is the first to prove ergodicity of a stochastic PDE in a hypoelliptic setting 
under conditions which compare favorably to those under which similar theo- 
rems are proven for finite dimensional stochastic differential equations. One of 
the keys to accomplishing this is a recent result from [MP06] on the regularity 
of the Malliavin matrix in this setting. 

One of the main technical contributions of the present work is to provide 
an infinitesimal replacement for Girsanov's theorem in the infinite dimensional 
nonadapted setting which the application of these ideas to the fully hypoelliptic 
setting seems to require. Another of the principal technical contributions is to 
observe that the strong Feller property is neither essential nor natural for the 
study of ergodicity in dissipative infinite-dimensional systems and to provide 
an alternative. We define instead a weaker asymptotic strong Feller property 
which is satisfied by the system under consideration and is sufficient to give 
ergodicity. In many dissipative systems, including the stochastic Navier-Stokes 
equations, only a finite number of modes are unstable. Conceivably, these 
systems are ergodic even if the noise is transmitted only to those unstable 
modes rather than to the whole system. The asymptotic strong Feller property 
captures this idea. It is sensitive to the regularization of the transition densities 
due to both probabilistic and dynamic mechanisms. 

This paper is organized as follows. In Section 2 the precise mathematical 
formulation of the problem and the main results for the stochastic Navier- 
Stokes equations are given. In Section 3 we define the asymptotic strong 
Feller property and prove in Theorem 3.16 that, together with an irreducibil- 
ity property it implies ergodicity of the system. We thus obtain the analog in 
our setting of the classical result often derived from theorems of Khasminskii 
and Doob which states that topological irreducibility, together with the strong 
Feller property, implies uniqueness of the invariant measure. The main tech- 
nical results are given in Section 4, where we show how to apply the abstract 
results to our problem. Although this section is written with the stochastic 
Navier-Stokes equations in mind, most of the corresponding results hold for a 
much wider class of stochastic PDEs with polynomial nonlinearities. 
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2. Setup and main results 

Consider the two-dimensional, incompressible Navier-Stokes equations on 
the torus T^ = [— 7r,7r]^ driven by a degenerate noise. Since the velocity and 
vorticity formulations are equivalent in this setting, we choose to use the vor- 
ticity equation as this simplifies the exposition. For u a divergence- free velocity 
field, we define the vorticity w hy w = V Au = d2Ui — diU2- Note that u can be 
recovered from w and the condition V -u = 0. With this notation the vorticity 
formulation for the stochastic Navier-Stokes equations is as follows: 

(2.1) dw = iyAwdt + B{ICw,w)dt + QdW{t) , 

where A is the Laplacian with periodic boundary conditions and B{u,w) = 
— {u-'V)w, the usual Navier-Stokes nonlinearity. The symbol QdW{t) denotes 
a Gaussian noise process which is white in time and whose spatial correlation 
structure will be described later. The operator /C is defined in Fourier space 
by {}Cw)k = —iwkk-^/\\k\\'^, where (fci,A;2)"'" = (^2,— ^i)- By Wk, we mean 
the scalar product of w with {2tt)~^ exp(i/c ■ x). It has the property that the 
divergence of ICw vanishes and that w = VA{ICw). Unless otherwise stated, we 
consider (2.1) as an equation inTi = Lq, the space of real- valued square-integra- 
ble functions on the torus with vanishing mean. Before we go on to describe 
the noise process QW, it is instructive to write down the two-dimensional 
Navier-Stokes equations (without noise) in Fourier space: 

,2 1 >r- /., A / 1 1 



(2.2) w, = -u\k\^w, - - ^ (^jU) [jm - M2 j^.-^^ • 

j+e=k ' ' '■'' 

From (2.2), we see clearly that any closed subspace of 7i spanned by Fourier 
modes corresponding to a subgroup of Z is invariant under the dynamics. In 
other words, if the initial condition has a certain type of periodicity, it will be 
retained by the solution for all times. 

In order to describe the noise QdW{t), we start by introducing a conve- 
nient way to index the Fourier basis of Tt. We write Z \ {(0, 0)} = Z_|_ U Z_, 
where 

Zl = {{ki,k2) e'Z?\k2 >0}u{(A;i,0) G Z^l/ei >0} , 

7?_ = {{ki,k2)ez^\ -kezl) , 
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(note that Z_[_ is essentially the upper half-plane) and set, for A; G Z \ {(0, 0)}, 

f sm(k ■ x) if /c G Zi, 
[ cos(A; ■ x) if A; G Zl. 

We also fix a set 

(2.4) Zo = {kn\n=l,...,m}cZ^\ {(0, 0)} , 

which encodes the geometry of the driving noise. The set Zq will correspond 
to the set of driven modes of equation (2.1). 

The process W{t) is an ?7i-dimensional Wiener process on a probabil- 
ity space {Q,,T,P). For definiteness, we choose Q, to be the Wiener space 
Co([0,oo),R"^), W the canonical process, and P the Wiener measure. We de- 
note expectations with respect to P by E and define J-t to be the a-algebra 
generated by the increments of W up to time t. We also denote by {e„} the 
canonical basis of R™. The linear map Q : R™ -^ 7i is given by Qsn = qnfk„, 
where the qn are some strictly positive numbers, and the wave numbers kn 
are given by the elements of Zq. With these definitions, QW is an "H-valued 
Wiener process. We also denote the average rate at which energy is injected 
into our system by Sq = trQQ* = ^„ g^. 

We assume that the set Zq is symmetric, i.e. that if A; G Zq, then —k G Zq. 
This is not a strong restriction and is made only to simplify the statements 
of our results. It also helps to avoid the possible confusion arising from the 
slightly nonstandard definition of the basis fk- This assumption always holds 
for example if the noise process QW is taken to be translation invariant. In 
fact. Theorem 2.1 below holds for nonsymmetric sets Zq if one replaces Zq in 
the theorem's conditions by its symmetric part. 

It is well-known [Fla94], [MR04] that (2.1) defines a stochastic flow on 7i. 
By a stochastic flow, we mean a family of continuous maps <I>i : ri x H ^ H 
such that wt = <I>t(VF,wo) is the solution to (2.1) with initial condition wq and 
noise W. Hence, its transition semigroup Vt given by Ptfiwo) = Ewgip{wt) is 
Feller. Here, ip denotes any bounded measurable function from "H to R and we 
use the notation E^^ for expectations with respect to solutions to (2.1) with 
initial condition wq. Recall that an invariant measure for (2.1) is a probability 
measure //* on Ti. such that V^fi-k = ^*, where V* is the semigroup on measures 
dual to Vt- While the existence of an invariant measure for (2.1) can be proved 
by "soft" techniques using the regularizing and dissipativity properties of the 
flow [Cru89], [Fla94], showing its uniqueness is a challenging problem that 
requires a detailed analysis of the nonlinearity. The importance of showing the 
uniqueness of ^+ is illustrated by the fact that it implies 

(2.5) lim — / (p{wt)dt= / (f{w)n^,{dw), 

T^OD 1 Jq Jj^ 
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for all bounded continuous functions ^ and ^^,-almost every initial condition 
1^0 G W. It thus gives some mathematical ground to the ergodic assumption 
usually made in the physics literature in a discusion of the qualitative behavior 
of (2.1). The main results of this article are summarized by the following 
theorem: 

Theorem 2.1. Let Zq satisfy the following two assumptions: 

Al. There exist at least two elements in Zq with different Euclidean norms. 

A2. Integer linear combinations of elements of Zq generate Z . 

Then, (2.1) has a unique invariant measure in 7i. 

Remark 2.2. As pointed out by J. Hanke, condition A2 above is equiva- 
lent to the easily verifiable condition that the greatest common divisor of the 
set { det{k,i) : k,i G Zq} is 1, where det{k,£) is the determinant of the 2x2 
matrix with columns k and l. 

The proof of Theorem 2.1 is given by combining Corollary 4.2 with Propo- 
sition 4.4 below. A partial converse of this ergodicity result is given by the 
following theorem, which is an immediate consequence of Proposition 4.4. 

Theorem 2.3. There are two qualitatively different ways in which the 
hypotheses of Theorem 2.1 can fail. In each case there is a unique invariant 
measure supported on 7i, the smallest closed linear subspace of 7i which is 
invariant under (2.1). 

• In the first case the elements of Zq are all collinear or of the same 
Euclidean length. Then Ti. is the finite- dimensional space spanned by 
{fk\k G Zq}, and the dynamics restricted to 7i is that of an Omstein- 
Uhlenbeck process. 

• In the second case let Q be the smallest subgroup of Z containing Zq. 
Then 7i is the space spanned by {fk | A; G ^ \ {(0,0)}}. Let ki, k2 be 
two generators for Q and define Vi = 27r ki/\ki\'^, then 7i is the space of 
functions that are periodic with respect to the translations vi and f 2 . 

Remark 2.4. That 7i constructed above is invariant is clear; that it is 
the smallest invariant subspace follows from the fact that the transition prob- 
abilities of (2.1) have a density with respect to the Lebesgue measure when 
projected onto any finite-dimensional subspace of 7i; see [MP06]. 

By Theorem 2.3 if the conditions of Theorem 2.1 are not satisfied then 
one of the modes with lowest wavenumber is in Ti . In fact either /(i,o) -L W or 
fn^i) -L Ti. On the other hand for sufficiently small values of u the low modes of 
(2.1) are expected to be linearly unstable [Fri95]. If this is the case, a solution 
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to (2.1) starting in TY"*" will not converge to TY and (2.1) is therefore expected 
to have several distinct invariant measures on TC. It is however known that the 
invariant measure is unique if the viscosity is sufficiently high; see [Mat99]. (At 
high viscosity, all modes are linearly stable. See [Mat03] for a more streamlined 
presentation.) 

Example 2.5. The set Zq = {(1,0), (-1, 0), (1, 1), (-1, -1)} satisfies the 
assumptions of Theorem 2.1. Therefore, (2.1) with noise given by 

QW{t, x) =Wi{t) sinxi + W2{t)cosxi + W3(t)sin(xi +X2) 
+W4(t)cos(xi +X2) , 

has a unique invariant measure in TC for every value of the viscosity i^ > 0. 

Example 2.6. Take Zq = {(1,0), (— 1, 0), (0, 1), (0, — 1)} whose elements 
are of length 1. Therefore, (2.1) with noise given by 

(2.6) QW{t,x) = Wi{t)smxi + W2{t)cosxi + Wsit) sin X2 + W^it) cos X2 , 

reduces to an Ornstein-Uhlenbeck process on the space spanned by sinxi, 
cosxi, sin 2:2, and cos2;2. 

Example 2.7. Take Zq = {(2,0), (-2, 0), (2, 2), (-2, -2)}, which corre- 
sponds to case 2 of Theorem 2.3 with Q generated by (0,2) and (2,0). In 
this case, H is the set of functions that are vr-periodic in both arguments. Via 
the change of variables x >-^ x/2, one can easily see from Theorem 2.1 that 
(2.1) then has a unique invariant measure on TC (but not necessarily on TC). 

3. An abstract ergodic result 

We start by proving an abstract ergodic result, which lays the foundations 
of the present work. Recall that a Markov transition semigroup Vt is said to 
be strong Feller at time t if Vt^p is continuous for every bounded measurable 
function ip. It is a well-known and much used fact that the strong Feller prop- 
erty, combined with some irreducibility of the transition probabilities implies 
the uniqueness of the invariant measure for Vt [DPZ96, Th. 4.2.1]. If Vt is 
generated by a diffusion with smooth coefficients on R" or a finite-dimensional 
manifold, Hormander's theorem [Hor67], [Hor85] provides us with an efficient 
(and sharp if the coefficients are analytic) criterion for the strong Feller prop- 
erty to hold. Unfortunately, no equivalent theorem exists if Vt is generated by 
a diffusion in an infinite-dimensional space, where the strong Feller property 
seems to be much "rarer" . If the covariance of the noise is nondegenerate (i.e. 
the diffusion is elliptic in some sense), the strong Feller property can often 
be recovered by means of the Bismut-Elworthy-Li formula [EL94]. The only 
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result to our knowledge that shows the strong Feller property for an infinite- 
dimensional diffusion where the covariance of the noise does not have a dense 
range is given in [EHOl], but it still requires the forcing to act in a nondegen- 
erate way on a subspace of finite codimension. 

3.1. Preliminary definitions. Let X he a Polish (i.e. complete, separable, 
metrizable) space. Recall that a pseudo-metric for A' is a continuous function 
d : X^ -^ R+ such that d{x, x) = and such that the triangle inequality is 
satisfied. We say that a pseudo-metric di is larger than ^2 if di{x, y) > d2{x, y) 
for ah {x,y) G X'^. 

Definition 3.1. Let {(i„}^g be an increasing sequence of (pseudo-)metrics 
on a Polish space X. If lim„^oo dn{x, y) = 1 for all x / y, then {d„} is a totally 
separating system of (pseudo-) metrics for X. 

Let us give a few representative examples. 

Example 3.2. Let {a„} be an increasing sequence in R such that 
lim„^oo On = oo. Then, {dn} is a totally separating system of (pseudo-) metrics 
for X in the following three cases. 

1. Let d be an arbitrary continuous metric on X and set dn{x,y) = 1 A 
and{x,y). 

2. Let X = Co(R) be the space of continuous functions on R vanishing at 
infinity and set dn{x,y) = 1 A sup^^r „„i an\x{s) — y{s)\. 

3. Let X = l"^ and set dn{x, y) = 1 A a„ Y.k=o \^k - Vkl"^- 

Given a pseudo-metric d, we define the following seminorm on the set of 
d-Lipschitz continuous functions from X to R: 

(3.1) ii^iu = -p MpiM . 

a,^y(zx d{x,y) 

This in turn defines a dual seminorm on the space of finite signed Borel mea- 
sures on X with vanishing integral by 



(3.2) III^IIU = sup / (p{x)iy{dx). 

Given ni and fi2, two positive finite Borel measures on X with equal mass, we 
also denote by C(/ii,/i2) the set of positive measures on X^ with marginals /ii 
and iJ,2 and we define 

(3.3) WfJ-i - IJ-2\\d = inf / d{x,y)fi{dx,dy) . 
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The following lemma is an easy consequence of the Monge-Kantorovich duality; 
see e.g. [Kan42], [Kan48], [AN87], and shows that in most cases these two 
natural notions of distance can be used interchangeably. 

Lemma 3.3. Let d be a continuous pseudo-metric on a Polish space X 
and let ni and fi2 be two positive measures on X with equal mass. Then, 

IImi ~ M2IU = IIImi ~ M2||U- 

Proof. This result is well-known if {X, d) is a separable metric space; see 
for example [Rac91] for a detailed discussion on many of its variants. If we 
define an equivalence relation on A" by 2; ~ y <^ d{x, y) = and set X^ = X j^, 
then d is well-defined on X^ and (^X^, d) is a separable metric space (although 
it may no longer be complete). When n : X ^ Xa hy 7r(x) = [x], the result 
follows from the Monge-Kantorovich duality in X^ and the fact that both sides 
of (3.3) do not change if the measures m are replaced by 7r*^j. D 

Recall that the total variation norm of a finite signed measure fi on X 
is given by ||^||tv = ^(m^('^) + /"~('^))) where ji = ^i^ — ijl~ \s the Jordan 
decomposition of ^. The next result is crucial to the approach taken in this 
paper. 

Lemma 3.4. Let {dn} he a hounded and increasing family of continuous 
pseudo-metrics on a Polish space X and define d{x,y) = lim„_>oo c^n(a;,y). 
Then, lim„__»oo ||/Ui — /U2||d„ = ||^i — /U2||d for any two positive measures fii and 
II2 with equal mass. 

Proof. The limit exists since the sequence is bounded and increasing 
by assumption, so let us denote this limit by L. It is clear from (3.3) that 
IImi ~ A*2||d > L, so it remains to show the converse bound. Let /i„ be a 
measure in Q(ni,H2) that realizes (3.3) for the distance d„. (Such a measure is 
shown to exist in [Rac91].) The sequence {/U„} is tight on X'^ since its marginals 
are constant, and so we can extract a weakly converging subsequence. Denote 
by /ioo the limiting measure. For m > n 

/ dn{x,y)firn{dx,dy)< / d„i{x,y)firn{dx,dy)<L. 

Since d„ is continuous, the weak convergence taking m ^ 00 implies that 



dn{x,y) Hoo{dx,dy) < L , Vn > . 

It follows from the dominated convergence theorem that J^2 d{x, y) fiooidx, dy) 
< L, which concludes the proof. D 

Corollary 3.5. Let X he a Polish space and let {dn} he a totally separat- 
ing system of pseudo-metrics for X . Then, ||^i — jU2||tv = lim„_»oo ||/Wi — ^2|U„ 
for any two positive measures fii and ^2 with equal mass on X . 
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Proof. It suffices to notice that 

ll/Ui -At2||TV = inf fi{{{x,y) :xy^y}) = \\iJ.i- H2\\d 
/^ee(/ii,/i2) 

with d{x, y) = 1 whenever x ^ y and then to apply Lemma 3.4. Observe 
that dn ^ dhj the definition of a totally separating system of pseudo-metrics 
and that Lemma 3.4 makes no assumptions on the continuity of the limiting 
pseudo-metric d. D 

3.2. Asymptotic strong Feller. Before we define the asymptotic strong 
Feller property, recall that: 

Definition 3.6. A Markov transition semigroup on a Polish space X is said 
to be strong Feller at time t if V^ip is continuous for every bounded measurable 
function ip : X ^ Ti. 

Note that if the transition probabilities Vtix, ■ ) are continuous in x in the 
total variation topology, then Vt is strong Feller at time t. 

Recall also that the support of a probability measure fi, denoted by 
supp(^), is the intersection of all closed sets of measure 1. A useful char- 
acterization of the support of a measure is given by 

Lemma 3.7. A point x £ supp(^) if and only if ^[U) > for every open 
set U containing x. 

It is well-known that if a Markov transition semigroup Vt is strong Feller 
and Hi and 112 are two distinct ergodic invariant measures for Vt (i.e. ni and 
H2 are mutually singular), then supp^i Pi supp/U2 = (p- (This can be seen 
e.g. by the same argument as in [DPZ96, Prop. 4.1.1].) In this section, we 
show that this property still holds if the strong Feller property is replaced by 
the following property, where we denote by Ux the collection of all open sets 
containing x. 

Definition 3.8. A Markov transition semigroup Vt on a Polish space X 
is called asymptotically strong Feller at x if there exists a totally separating 
system of pseudo-metrics {dn} for X and a sequence tn > such that 

(3.4) inf lim sup sup ||Pt„(x, • ) - VtAv, • )IU„ = , 

It is called asymptotically strong Feller if this property holds at every x £ X. 

Remark 3.9. If ^6(2;, 7) denotes the open ball of radius 7 centered at x 
in some metric defining the topology of X, then it is immediate that (3.4) is 
equivalent to 

lim lim sup sup \\Vt„ix,-) - Vt„{y,-)\\d„ = . 

7-^0 n^oo yeB{x,-t) 
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Remark 3.10. Notice that the definition of the asymptotic strong Feller 
property allows for the possibility that tn = t for all n. In this case, the 
transition probabilities Vtix,-) are continuous in the total variation topology 
and thus Vg is strong Feller at times s > t. Conversely, strong Feller Markov 
semigroups on Polish spaces are asymptotically strong Feller. To see this first 
observe that if P and Q are two Markov operators over the same Polish space 
that are strong Feller, then the product PQ is a Markov operator whose transi- 
tion probabilities are continuous in the total variation distance [DM83], [Sei02]. 
Hence, if Vt is strong Feller for some t > 0, then 7^2* = PtPt is continuous in 
the total variation distance, which implies that the semigroup Vt is asymptot- 
ically strong Feller. We would like to thank B. Goldys for pointing this fact 
out to us. 

One other way of seeing the connection to the strong Feller property is to 
recall that a standard criterion for Vt to be strong Feller is given by [DPZ96, 
Lem. 7.1.5]: 

Proposition 3.11. A semigroup Vt on a Hilhert space 7i is strong Feller 
if, for all if : 7i ^ H with Hv^Hoo = sup^gT^ Iv^l^^)! ^^t? ||V(/3||oo finite one has 

(3.5) |VPt^(x)| < C{\\x\\)M^ , 

where C : R+ ^ H is a fixed nondecreasing function. 

The following lemma provides a similar criterion for the asymptotic strong 
Feller property: 

Proposition 3.12. Let tn and 5n be two positive sequences with {tn} 
nondecreasing and {6n} converging to zero. A semigroup Vt on a Hilhert space 
Ti is asymptotically strong Feller if, for all ip : Ti. ^ H with ||9?||oo and \\Vip\\oo 
finite, 

(3.6) NVtMx)\ < C{\\x\\){\\ip\\^ + 5n\\Vip\U) 
for all n, where C : R+ -^ Ti is a fixed nondecreasing function. 

Proof. For e > 0, we define on 7i the distance 

d£{vji,W2) = 1 A e^-^ll^i - W2\\, 

and we denote by || • ||e the corresponding seminorms on functions and on 
measures given by (3.1) and (3.2). It is clear that if 5n is a decreasing sequence 
converging to 0, {ds^} is a totally separating system of metrics for 7i. 
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It follows immediately from (3.6) that for every Frechet differentiable func- 
tion ip from "H to R with \\'p\\e < 1, 
(3.7) 

if{w){VtAwudw)-VtSw2,dw)) < \\wi - W2\\C{\\wi\\V \\w2\\)(l + —) . 
n V e / 

Now take a Lipschitz continuous function ip with \\f\\e < 1- By apply- 
ing to f the semigroup at time 1/m corresponding to a linear strong Feller 
diffusion in 7i, one obtains ([Cer99], [DPZ96]) a sequence prn of Frechet differ- 
entiable approximations ipm with ||93m||e ^ 1 and such that (fm -^ f pointwise. 
Therefore, by the dominated convergence theorem, (3.7) holds for Lipschitz 
continuous functions (p and so 

WVtAwi, • ) - rtAw2, ■)\\e< \\WI - W2\\C{\\wi\\ V ||u;2||) (l + ^) • 

Choosing e = a„ = \/3^, we obtain 

\\VtAwi,-)-'PtAw2,-)\\a^ < ||ti'i-ti'2||C(||u'i|| V||u;2||)(l + a„) , 
which in turn implies that Vt is asymptotically strong Feller since a„ ^ 0. D 

Example 3.13. Consider the SDE 

dx = —xdt + dW{t) , dy = —y dt . 

Then, the corresponding Markov semigroup Vt on R is not strong Feller, 
but it is asymptotically strong Feller. To see that Vt is not strong Feller, 
let ip{x,y) = sgn(y) and observe that Vtp> = p for all t G [0, oo). Since p 
is bounded but not continuous, the system is not strong Feller. To see that 
the system is asymptotically strong Feller observe that for any differentiable 
9? : R ^ R and any direction ^ G R with ||^|| = 1, 

\{Wt(p){xo,yo) ■ ^\ = \E(^^^^y^){Vip){xt,yt) ■ {ut,vt)\ 

<||V(^||ooE|(nt,wt)| < \\Vip\\^e~\ 

where {ut,vt) is the linearized flow starting from ^. In other words {uq, vq) = ^, 
du = —udt, and dv = —vdt. This is a particularly simple example because the 
flow is globally contractive. 

Example 3.14. Now consider the SDE 

dx = {x - x^) dt + dW{t) , dy = -ydt . 

Again the function (p{x, y) = sgn(y) is invariant under Vt implying that the 
system is not strong Feller. It is however not globally contractive. As in the 
previous example, let ^ = (^1,^2) G R^ and ||.^|| = 1 and now let {ut,vt) 
denote the linearizion of this equation with (uq,vo) = C- Let "Pf denote the 
Markov transition semigroup of the xt process. It is a classical fact that for 
such a uniformly elliptic diffusion with a unique invariant measure one has 
\dxVt<pix,y)\ < C(|x|)||99||oo for some nondecreasing function C and all i > 1. 
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Hence differentiating witli respect to botli initial conditions produces 

<C(|x|)||(^||oo+E|?;t|||V(/?||oo 
<(C(|x|) + l)(||(^|U + e-*||Vv.||oo) 

for t > 1 which imphes that the system is asymptotically strong Feller. 

Example 3.15. In infinite dimensions, even a seemingly nondegenerate dif- 
fusion can suffer from a similar problem. Consider the following infinite di- 
mensional Ornstein-Uhlenbeck process u{x,t) = ^ii(A;, t) exp(iA;x) written in 
terms of its complex Fourier coefficients. We take x G T = [— tt, tt], k £ Z and 

(3.8) du{k,t) = -{1 + \k\'^)u{k,t) dt + exp{-\kf) dpkit) , 

where the Pk are independent standard complex Brownian motions. The 
Markov transition densities Vtix, ■) and Vt{y, ■) are singular for all finite times 
if X — y is not sufficiently smooth. This implies that the diffusion (3.8) in 
Tt = L^([— TTjVr]) is not strong Feller since by Lemma 7.2.1 of [DPZ96] the 
strong Feller property is equivalent to Vt{y, •) being equivalent to Vt{x, •) for 
all X and y. Another equivalent characterization of the strong Feller property is 
that the image(S'i) C image((5i) where St is the linear semigroup generated by 
the deterministic part for the equation defined by {Stu){k) = e~^^^'^' '^u{k,0) 
and Qt = L SrGS*dr where G is the covariance operator of the noise defined 
by {Gu){k) = exp{—2\k\^)u{k) . This captures the fact that the mean, con- 
trolled by St, is moving towards zero too slowly relative to the decay of the 
noise's covariance structure. However, one can easily check that the example 
is asymptotically strong Feller since the entire flow is pathwise contractive as 
in the first example. 

The classical strong Feller property captures well the smoothing due to 
the random effects. When combined with irreducibility in the same topology, 
it implies that the transition densities starting from different points are mu- 
tually absolutely continuous. As the examples show, this is often not true in 
infinite dimensions. We see that the asymptotic strong Feller property better 
incorporates the smoothing due to the pathwise contraction of the dynamics. 
Comparing Proposition 3.11 with Proposition 3.12, one sees that the second 
term in Proposition 3.12 allows one to capture the progressive smoothing in 
time from the pathwise dynamics. This becomes even clearer when one ex- 
amines the proofs of Proposition 4.3 and Proposition 4.11 later in the text. 
There one sees that the first term comes from shifting a derivative from the 
test function to the Wiener measure and the second is controlled using in an 
essential way the contraction due to the spatial Laplacian. 
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The usefulness of the asymptotic strong Feller property is seen in the 
following theorem and its accompanying corollary which are the main results 
of this section. 

Theorem 3.16. Let Vt be a Markov semigroup on a Polish space X and 
let fi and v he two distinct ergodic invariant probability measures for Vt- IfVt 
is asymptotically strong Feller at x, then x supp^ n suppz^. 

Proof. By Corollary 3.5, the proof of this result is a simple rewriting of 
the proof of the corresponding result for strong Feller semigroups. 

For every measurable set A, every t > 0, and every pseudo-metric d on X 
with d < 1, the triangle inequality for || • H^^ implies 



(3.9) 11// - z/|U < 1 - mm{fiiA),u{A)}[l - max \\Vtiz, •) - Vtiy, ■)\\d 

To see this, set a = inm{fi{A),i'{A)}. If a = there is nothing to prove 
so assume a > 0. Clearly there exist probability measures P, /Z, i^a, and 
fj,A such that i^A^A) = HAiA) = 1 and such that n = {1 — a)fl + ajiA and 
V = [1 — a)i> + avA- From the invariance of the measures (i and f and the 
triangle inequality this implies 

IIa* - T^Wd = WVtfJ, - Vfi^Wd < (1 - a)\\Vtfl - Vti>\\d + a\\VtfJ-A - Vt^AWd 

< (1 - a) + Q / / \\Vt{z, •) - Vt{y, ■)\\diiA{dz)uA{dy) 
JaJa 



< 1-a 1- max WVtiz,-)- Vtiy, ■)\ 
\ y,zeA 

Continuing with the proof of the corollary, we see that, by the definition of 
the asymptotic strong Feller property, there exist constants iV > 0, a sequence 
of totally separating pseudo-metrics {dn}, and an open set U containing x such 
that \\Vt„{z, •) — Vt^{y, •)\\d„ < 1/2 for every n > N and every y,z £ U. (Note 
that by the definition of totally separating pseudo-metrics dn < 1.) 

Assume by contradiction that x £ supp/i n suppz/ and therefore that 
a = iam{fi{U), v{U)) > 0. Taking A = U, d = dn, and t = t„ in (3.9), we then 
get Wfi — i^\\d„ < 1 — ^ for every n > N, and therefore \\fi — z^||tv < 1 — f by 
Corollary 3.5, thus leading to a contradiction. D 

As an immediate corollary, we have 

Corollary 3.17. IfVt is an asymptotically strong Feller Markov semi- 
group and there exists a point x such that x £ supp/i for every invariant 
probability measure fi ofPt, then there exists at most one invariant probability 
measure for Vt ■ 
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4. Applications to the stochastic 2D Navier-Stokes equations 

To state the general ergodic result for the two-dimensional Navier-Stokes 
equations, we begin by looking at the algebraic structure of the Navier-Stokes 
nonlinear ity in Fourier space. 

Remember that Zq as given in (2.4) denotes the set of forced Fourier 
modes for (2.1). In view of Equation 2.2, it is natural to consider the set ^oc 
defined as the smallest subset of Z containing Z^^ and satisfying that for every 
£, j G Joo such that (^^, j) / and |j| / \i\, one has j + ^ G ^oo (see [EMOl]). 
Denote by li. the closed subspace of H. spanned by the Fourier basis vectors 
corresponding to elements of Z^a- Then, H, is invariant under the flow defined 
by (2.1). 

Since we would like to make use of the existing results, we recall the 
sequence of subsets Z^ of T? defined recursively in [MP06] by 

J G ^0,^ G ^„_iwith U^.j) / 0,|j| / l^ll , 



Zn = \i^3 

as well as Zoo = U^i ^n- The two sets Z^o and Z^o are the same even though 
from the definitions we only see Z^o C Z^o- The other inclusion follows from 
the characterization of Z^o given in Proposition 4.4 below. 

The following theorem is the principal result of this article. 

Theorem 4.1. The transition semigroup onTC generated by the solutions 
to (2.1) is asymptotically strong Feller. 

An almost immediate corollary of Theorem 4.1 is 

Corollary 4.2. There exists exactly one invariant probability measure 
for (2.1) restricted to 7i. 

Proof of Corollary 4.2. The existence of an invariant probability measure 
H for (2.1) is a standard resuh [Fla94], [DPZ96], [CK97]. By Corollary 3.17 
it suffices to show that the support of every invariant measure contains the 
element 0. Applying Ito's formula to ||^«|p yields for every invariant measure 
/i the a priori bound 

I l|2 /J \ ^ ^"^0 

\w\\ fi{aw) < . 

H ^ 

(See [EMSOl, Lemma B.l].) Therefore, denoting by Bi^p) the ball of radius p 
centered at 0, we have C such that \x{B((J^ > ^ for every invariant measure /x. 
On the other hand, [EMOl, Lemma 3.1] shows that, for every 7 > there exists 
a time Ty such that 

inf Vt,{w,B{j)) >0. 

weB{C) 
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(Note, though [EMOl, Lemma 3.1] was about Galerkin approximations, inspec- 
tion of the proof reveals that it holds equally for the full solution.) Therefore, 
fi{B{^)) > for every 7 > and every invariant measure ^, which implies that 
G supp(;u) by Lemma 3.7. D 

The crucial ingredient in the proof of Theorem 4.1 is the following result: 

Proposition 4.3. For every ry > 0, there exist constants C,5 > such 
that for every Frechet dijjerentiable function cp from 7i toTi one has the bound 

(4.1) \\VVnV{w)\\ < CeMv\\wf){M\oo + ||V(/.||ooe-^") , 

for every w £ 7i and n G N. 

The proof of Proposition 4.3 is the content of Section 4.6 below. Theo- 
rem 4.1 then follows from this proposition and from Proposition 3.12 with the 
choices tn = n and dn = e~ . Before we turn to the proof of Proposition 4.3, 
we characterize Z^o and give an informal introduction to Malliavin calculus 
adapted to our framework, followed by a brief discussion on how it relates to 
the strong Feller property. 

4.1. The structure of Zoo- In this section, we give a complete characteri- 
zation of the set Zoo- We start by defining (Zq) as the subset of Z^ \ {(0, 0)} 
generated by integer linear combinations of elements of Zq. With this notation, 
we have 

Proposition 4.4. If there exist 01,02 G Zq such that \ai\ 7^ \a2\ and such 
that ai and 02 are not coUinear, then Z^o = {Zq). Otherwise, Z^o = Zq. In 
either case, one always has that Z^o = Z^ 



-'CO- 



This also allows us to characterize the main case of interest: 

Corollary 4.5. One has Z^o = Z^ \ {(0, 0)} if and only if the following 
holds: 

1. Integer linear combinations of elements of Zq generate Z . 

2. There exist at least two elements in Zq with nonequal Euclidean norm- 
Proof of Proposition 4.4. It is clear from the definitions that if the 

elements of Zq are all coUinear or of the same Euclidean length, one has Z^o = 
Zq = Zoo- In the rest of the proof, we assume that there exist two elements 
a\ and 02 of Zq that are neither collinear nor of the same length and we 
show that one has Z^o = {Zq). Since it follows from the definitions that 
Zoo C Zoo C {Zq) , this shows that Zoo = ^oo - 
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Note that the set Z^o consists exactly of those points in Z that can be 
reached by a walk starting from the origin with steps drawn in Zq and which 
does not contain any of the following "forbidden steps" : 



Definition 4.6. A step with increment i ^ Zq starting from j G 
forbidden if either \j\ = \i\ or j and i are collinear. 



IS 



Our first aim is to show that there exists R > such that Z^o contains 
every element of {Zq) with Euclidean norm larger than R. In order to achieve 
this, we start with a few very simple observations. 

Lemma 4.7. For every Rq > 0, there exists Ri > such that every j G 
(Zq) with \j\ < Rq can be reached from the origin by a path with steps in Zq 
{some steps may be forbidden) which never exits the ball of radius Ri. 

Lemma 4.8. There exists L > such that the set Z^q contains all ele- 
ments of the form niOi + n2a2 with rii and n2 m Z \ [—L, L]. 

Proof. We may assume without loss of generality that |ai| > \a2\ and 
that (ai,a2) > 0. Choose L such that L (01,02) > |oiP- By the symmetry 
of Zq, we can replace (01,02) by (—01,-02), so that we can assume without 
loss of generality that 77-2 > 0. We then take first one step in the direction oi 
starting from the origin, followed by 77-2 steps in the direction 02. Note that 
the assumptions we made on oi, 02, and n2 ensure that none of these steps is 
forbidden. From there, the condition n2 > L ensures that we can take as many 
steps as we want into either the direction oi or the direction — oi without any 
of them being forbidden. D 

Denote by Z the set of elements of the form nioi + 1x202 considered in 
Lemma 4.8. It is clear that there exists Rq > such that every element in 
(Zq) is at distance less than Rq of an element of Z. Given this value Rq, we 
now fix i?i as given from Lemma 4.7. Let us define the set 

A = Z^n {{aj I Q G R, i G Zo} U {/c I 3i G Zq with|j| = \k\}) , 

which has the property that there is no forbidden step starting from Z \ A. 
Define furthermore 

B = {je{Zo) I mi\k-j\>Ri}. 

keA 

By Lemma 4.7 and the definition of B, every element of B can be reached 
by a path from Z containing no forbidden steps, therefore B C Z^o- On the 
other hand, it is easy to see that there exists R > such that for every element 
of J G {Zq) \B with IjI > R, there exists an element a{j) G Zq and an element 
k{j) G B such that j can be reached from k{j) with a finite number of steps in 
the direction a(j). Furthermore, if R is chosen sufficiently large, none of these 
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Figure 1: Construction from the proof of Proposition 4.4. 



steps crosses A, and therefore none of them is forbidden. We have thus shown 
that there exists R > such that Z^o contains {j G {Zq) | |jp > -R}. 

In order to help in visualizing this construction, Figure 1 shows the typical 
shapes of the sets A (dashed hues) and B (gray area), as weh as possible choices 
of a{j) and k{j), given j. (The black dots on the intersections of the circles 
and the hues making up A depict the elements oi Zq.) 

We can (and will from now on) assume that R is an integer. The last step 
in the proof of Proposition 4.4 is 

Lemma 4.9. Assume that there exists an integer R > 1 such that Z^ 



contains {j G {Zq) 
R-l}. 



\j\ > R}. Then Z^o also contains {j G {Zq) \\j\ > 



Proof. Assume that the set {j £ {Zq) | |jp = R — 1} is nonempty and 
choose an element j from this set. Since Zq contains at least two elements 
that are not collinear, we can choose k £ Zq such that k is not collinear to j. 
Since Zq is closed under the operation k i-^ —k, we can assume that (j, k) > 0. 
Consequently, one has \j + k\'^ > R, and so j + fc G Zc^ by assumption. The 



same argument shows that |j + k\ > \k\ + 1, so the step 
j + k is not forbidden and therefore k G Z^q. 



-k starting from 
D 



This shows that Zoo 
sition 4.4. 



{Zq) and therefore completes the proof of Propo- 



D 
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4.2. Malliavin calculus and the Navier-Stokes equations. In this section, 
we give a brief introduction to some elements of Malliavin calculus applied to 
equation (2.1) to help orient the reader and fix notation. We refer to [MP06] 
for a longer introduction in the setting of equation (2.1) and to [Nua95], [Bel87] 
for a more general introduction. 

Recall from Section 2, that ^t'- Ci[0,t]; R™") x 7i ^ 7i was the map such 
that wt = ^t(W,wo) for initial condition wq and noise realization W. Given 
a w G L^Qj,(R+, R™), the Malliavin derivative of the 7^-valued random variable 
Wt in the direction v, denoted T)'"wt, is defined by 

^t{W + eV,wo)-MW,wo) 
D Wt = hm 



-0 e 

where the limit holds almost surely with respect to the Wiener measure and 
where we set V{t) = Jq v{s) ds. Note that we allow v to be random and possibly 
nonadapted to the filtration generated by the increments of W. 

Defining the symmetrized nonlinearity B(w,v) = B{ICw,v) + B{ICv,w), 
we use the notation Jg^t with s < t for the derivative flow between times s 
and t, i.e. for every ^ £ 7i, Js^tC is the solution of 

(4.2) dtJs4 = i^AJs4 + Biwt,Js4), t>s, Js,s^ = ^. 



m\ 



Note that we have the important cocycle property J^f = Jr,tJs,r for r £ [s,t]. 

Observe that 1)'"wt = A^^tv where the random operator Ag^t '■ ^{[sj Aj'^ 
— > W is given by 

As,tv = / Jr,tQv{r) dr . 



To summarize, Jo^ti is the effect on wt of an infinitesimal perturbation of the 
initial condition in the direction ^ and AQ^tv is the effect on wt of an infinitesimal 
perturbation of the Wiener process in the direction of V{s) = L v{r) dr. 

Two fundamental facts we will use from Malliavin calculus are embodied 
in the following equalities. The first amounts to the chain rule, the second is 
integration by parts. For a smooth function (p : 7i ^ Ti and a (sufficiently 
regular) process v, 

(4.3) B{{Vip){wt),'D-wt) = B(v\^{wt))) ='E.(ip{wt) j {v{s),dWs 

The stochastic integral appearing in this expression is an Ito integral if the 
process v is adapted to the filtration Tt generated by the increments of W and 
a Skorokhod integral otherwise. 

We also need the adjoint A*^^ -.H -^ L^([s, i], R™") defined by the duality 
relation l^A*^ ^^, f ) = (^, Ag^tv), where the first scalar product is in L^([s, t], R™") 
and the second one is in Ti. Note that one has (yl*^^)(r) = Q*J*^^, where J*^ 
is the adjoint in Ti oi Jrt- 
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One of the fundamental objects in the study of hypoelhptic diffusions is 
the Malhavin matrix Mg^t = ^s t^*s t- ^ ghmpse of its importance can be seen 
from the following. For ^ G H, 

m „t 



'" /■t 



Hence the quadratic form (Mo,f ^, is zero for a direction ^ only if no variation 
whatsoever in the Wiener process at times s <t could cause a variation in wt 
with a nonzero component in the direction ^. 

We also recall that the second derivative Kg^t of the flow is the bilinear 
map solving 

dtKsM^ O = uAKsA^, O + B{wt, Ks,t{^, e')) + B{Js,t^', Js,tO , 

It follows from the variation-of-constants formula that Kg^ti^,^') is given by 
(4.4) Ks,t{^, CO = / Jr,tB{Js,re, JsA) dr ■ 

J s 

4.3. Motivating discussion. It is instructive to proceed formally pretending 
that Mq^i is invertible as an operator on 7i. This is probably not true for the 
problem considered here and we will certainly not attempt to prove it in this 
article, but the proof presented in Section 4.6 is a modification of the argument 
in the invertible case and hence it is instructive to start there. 

Set C( = Jo,tC; iiow C( can be interpreted as the perturbation of wt caused 
by a perturbation ^ in the initial condition of wt- Our goal is to find an 
infinitesimal variation in the Wiener path W over the interval [0, t] which 
produces the same perturbation at time t as the shift in the initial condition. 
We want to choose the variation which will change the value of the density the 
least. In other words, we choose the path with the least action with respect to 
the metric induced by the inverse of the Malliavin matrix. The least squares 
solution to this variational problem is easily seen to be, at least formally, 
V = Al^M^J£,t where v G L2([0, t], R™). Observe that Ti^wt = A^^tv = Jo,tC- 
Considering the derivative with respect to the initial condition w of the Markov 
semigroup Vt acting on a smooth function (p, we obtain 



(4.5) {VVMw),i) = E^((Vv9)(u;i)Jo,iO = E^((Vv9)(u;i)D"u)t) 



^^{(^{wt) I v{s)dWs\ < llv'llooE^ 



t 
v[s)dWs 





where the penultimate estimate follows from the integration by parts formula 
(4.3). Since the last term in the chain of implications holds for functions which 
are simply bounded and measurable, the estimate extends by approximation to 
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that class of if. Furthermore since the constant E«,|L t;(s)(iWs| is independent 
of (/3, if one can show it is finite and bounded independently oi ^ £ 7i with 
ll^ll = 1, we have proved that ||V'Pi(/3|| is bounded and thus that Vt is strong 
Feller in the topology of TC. Ergodicity then follows from this statement by 
means of Corollary 3.17. In particular, the estimate in (4.1) would hold. 

In slightly different language, since v is the infinitesimal shift in the 
Wiener path equivalent to the infinitesimal variation in the initial condition ^, 
we have, via the Cameron-Martin theorem, the infinitesimal change in the 
Radon-Nikodym derivative of the "shifted" measure with respect to the origi- 
nal Wiener measure. This is not trivial since in order to compute the shift v, 
one uses information on {ws}se[o,t]: so that it is in general not adapted to 
the Wiener process Ws- This nonadaptedness can be overcome as Section 4.8 
demonstrates. However the assumption in the above calculation that Mo,t 
is invertible is more serious. We will overcome this by using the ideas and 
understanding which begin in [Mat98], [Mat99], [EMSOl], [KSOO], [BKLOl]. 

The difficulty in inverting Mq^i partly lies in our incomplete understanding 
of the natural space in which (2.1) lives. The knowledge needed to identify on 
what domain Mo,t can be inverted seems equivalent to identifying the correct 
reference measure against which to write the transition densities. By "refer- 
ence measure," we mean a replacement for the role of Lebesgue measure from 
finite dimensional diffusion theory. This is a very difficult proposition. An al- 
ternative was given in the papers [Mat98], [Mat99], [KSOO], [EMSOl], [BKLOl], 
[Mat02b], [BKL02], [Hai02], [MY02]. The idea was to use the pathwise con- 
tractive properties of the flow at small scales due to the presence of the spatial 
Laplacian. Roughly speaking, the system has finitely many unstable directions 
and infinitely many stable directions. One can then use the noise to steer the 
unstable directions together and let the dynamics cause the stable directions 
to contract. This requires the small scales to be enslaved to the large scales 
in some sense. A stochastic version of such a determining modes statement 
(cf. [FP67]) was developed in [Mat98]. Such an approach to prove ergodicity 
requires looking at the entire future to -l-oo (or equivalently the entire past) 
as the stable dynamics only brings solutions together asymptotically. In the 
first papers in the continuous time setting ([EMSOl], [Mat02b], [BKL02]), Gir- 
sanov's theorem was used to bring the unstable directions together completely; 
[Hai02] demonstrated the effectiveness of only steering all of the modes together 
asymptotically. Since all of these techniques used Girsanov's theorem, they re- 
quired that all of the unstable directions be directly forced. This is a type of 
partial ellipticity assumption, which we will refer to as "effective ellipticity." 
The main achievement of this text is to remove this restriction. We also make 
another innovation which simplifies the argument considerably. We work in- 
finitesimally, employing the linearization of the solution rather than looking at 
solutions starting from two different starting points. 
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4.4. Preliminary calculations and discussion. Throughout this and the 
following sections we fix once and for all the initial condition wq & 7i for (2.1) 
and denote by wt the stochastic process solving (2.1) with initial condition 
wq. By E we mean the expectation starting from this initial condition unless 
otherwise indicated. Recall also the notation £q = tiQQ* = XlkfcP- ^^^^ 
following lemma provides the auxiliary estimates which will be used to control 
various terms during the proof of Proposition 4.3. 

Lemma 4.10. The solution of the 2D Navier-Stokes equations in the vor- 
ticity formulation (2.1) satisfies the following hounds: 

1. There exist positive constants C and r/o, depending only on Q and v, such 
that 

(4.6) Eexplrysupl ||tt)(|| + i^ \\wr\\i dr — £Q{t — s) 

^ t>s ^ Js 

<Cexp(r?e-''"||u;of) , 

for every s > and for every f] < rjo- Here and in the sequel, the notation 

\\w\\i = \\Vw\\ is used. 

2. There exist constants ryi,a, 7 > 0, depending only on £0 and v, such that 

N 

(4.7) Eexpl ry >^ ||tt)„|| — 7A''j < exp(ary||ti;o|| ) , 

n=0 

holds for every N > 0, every rj < rji, and every initial condition wq & TC. 

3. For every rj > 0, there exists a constant C = C{£Q,u,rj) > such that 
the Jacobian Jo,t satisfies almost surely 

(4.8) ||Jo,J < exp(^r/ / \\ws\\lds + Ct^ , 

for every t > 0. 

4. For every rj > and every p > 0, there exists C = C{£oi ^^ ViP) > such 
that the Hessian satisfies 

E\\K,4P < C exp{7]\\wof) , 

for every s > and every t £ (s, s + 1). 

The proof of Lemma 4.10 is postponed to Appendix A. 

We now show how to modify the discussion in Section 4.3 to make use of 
the pathwise contractivity on small scales to remove the need for the Malliavin 
covariance matrix to be invertible on all of Ti. 
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The point is that since the Malhavin matrix is not invertible, we are not 
able to construct a v £ L^([0,T],R™) for a fixed value of T that produces the 
same infinitesimal shift in the solution as an (arbitrary but fixed) perturbation 
^ in the initial condition. Instead, we will construct a v G L^([0,oo),R'") such 
that an infinitesimal shift of the noise in the direction v produces asymptot- 
ically the same effect as an infinitesimal perturbation in the direction ^. In 
other words, one has || Jo,tC ~ ^o,t^o,i|| ^ as t ^ oo, where fo,t denotes the 
restriction of v to the interval [0, t]. 

Set pt = Jo,tC ~ ^o,tVo,t, the residual error for the infinitesimal variation 
in the Wiener path W given by v. Then we have from (4.3) the approximate 
integration by parts formula: 



(4.9) {vrMw),o = ^w[{v{vim)),0) =^-'[i^v)(^t)Jo,tC 

= E^ (^ (V93) {wt)Ao,tVo,tj + E^ ( (Vv?) {wt)pt] 



■■^w[^iwt) f vis)dWis)^ +B^{{V^){wt)pt) 
ft 



< ll^'llooE^ 



v{s)dW{s) + ||V(^||ooE^||pj| 




This formula should be compared with (4.5). Again if the process v is not 
adapted to the filtration generated by the increments of the Wiener process 
W{s), the integral must be taken to be a Skorokhod integral; otherwise Ito 
integration can be used. Note that the residual error satisfies the equation 

(4.10) dtpt = i^Apt + B{wt, pt) - Qv{t) , po = C, 

which can be interpreted as a control problem, where v is the control and \\pt\\ 
is the quantity that one wants to drive to 0. 

If we can find a f so that pt ^ as t ^ co and E|L°° v{s) dW{s)\ < 00 
then (4.9) and Proposition 3.12 would imply that wt is asymptotically strong 
Feller. A natural way to accomplish this would be to take v{t) = Q~^B{wt, pt), 
so that dtpt = i^Apt and hence pt ^ as t ^ 00. However for this to make 
sense it would require that B{wt,pt) takes values in the range of Q. If the 
number of Brownian motions m is finite this is impossible. Even if m = 00, this 
is still a delicate requirement which severely limits the range of applicability 
of the results obtained (see [FM95], [Fer97], [MS05]). 

To overcome these difficulties, one needs to better incorporate the path- 
wise smoothing which the dynamics possesses at small scales. Though our 
ultimate goal is to prove Theorem 4.1, which covers (2.1) in a fundamentally 
hypoelliptic setting, we begin with what might be called the "essentially ellip- 
tic" setting. This allows us to outline the ideas in a simpler setting. 
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4.5. Essentially elliptic setting. To help to clarify the techniques used in 
the sections which follow and to demonstrate their applications, we sketch the 
proof of the following proposition which captures the main results of the earlier 
works on ergodicity, translated into the framework of the present paper. 

Proposition 4.11. Let Vt denote the semigroup generated by the solu- 
tions to (2.1) on TC. There exists an iV* = N^,{£q,u) such that if Zq contains 
{fceZ ,0<|A;|<A'^*}, then for any rj > there exist positive constants c 
and 7 so that 



\VVMw)\ < cexp (7]\\w\\'') ( ||99||oo + e-^iV(^|| 



This result translates the ideas in [EMSOl], [Mat02b], [Hai02] to our 
present setting. (See also [Mat03] for more discussion.) The result does differ 
from the previous analysis in that it proceeds infinitesimally. However, both 
approaches lead to proving the system has a unique ergodic invariant measure. 

The condition on the range of Q can be understood as a type of "effec- 
tive ellipticity." We will see that the dynamics is contractive for directions 
orthogonal to the range of Q. Hence if the noise smooths in these directions, 
the dynamics will smooth in the other directions. What directions are con- 
tracting depends fundamentally on a scale set by the balance between Sq and 
u (see [EMSOl, Mat03]). Proposition 4.3 holds given a minimal nondegener- 
acy condition independent of the viscosity ly, while Proposition 4.11 requires a 
nondegeneracy condition which depends on i'. 

Proof of Proposition 4.11. Let tt/j be the orthogonal projection onto the 
span of {fk : \k\ > N} and tt^ = 1 — tt/i. We will fix N presently; however, we 
will proceed assuming 7ii = TTiTC C Range (Q) and that Qe = ttiQ is invertible 
on Tif. By (4.10) we therefore have full control on the evolution of nept by 
choosing v appropriately. This allows for an "adapted" approach which does 
not require the control v to use information about the future increments of the 
noise process W. 

Our approach is first to define a process (t with the property that vr^Ct 
is after a finite time and iThCt evolves according to the linearized evolution, 
and then choose v such that pt = Ct- Since vr^Ct = after some time and the 
linearized evolution contracts the high modes exponentially, we readily obtain 
the required bounds on moments of pt- One can in fact pick any dynamics 
which are convenient for the modes which are directly forced. In the case 
when all of the modes are forced, the choice Ct = (1 — t/T)Jo^tC for t £ [0,T] 
produces the well-known Bismut-Elworthy-Li formula [EL94]. However, this 
formula cannot be applied in the present setting as all of the modes are not 
necessarily forced. 
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(4.11 



For ^ G H with ||^|| = 1, define (t by 
1 T^eCt 



dtCt 



+ uATThCt + 'n:hB{wt,Ct) , Co = C • 



2||vr,Ct|l 

With the convention that 0/0 = 0, set Cf = T^hCt and (f = Tr^Ct- We define the 
infinitesimal perturbation v by 

(4.12) v{t)=Qj'Ft, Ft = l-^ + uACf + 7r,B{wt,Ct)- 

^ IKt II 

Because Ft G H^, Q~^ Ft is well defined. It is clear from (4.10) and (4.12) 
that pt and Q satisfy the same equation, so that indeed pt = Ct- Since C,^ 



satisfies dt\\C-, 



t\\2 



-||Ci||, one has ||Ct || < ||Coll < l^ll = 1- Furthermore, for 



any initial condition wq and any ^ with ||^|| = 1, one has \\Q\\ = for t > 2. 
By calculations similar to those in Appendix A, there exists a constant C so 
that for any r/ > 



dtK 



h\\2 



< 



vN 



2 ^ II ||2\ n^h\\2 , ^ II W2,\\/- 

7Z7> - mwt\\i)\\Ct II +-lktllillC 



vrj^ 



V 



2||A^||2 



Hence, 

||Cf|P<||Co'lPexp 

+Cexp 



vN' 



C 
vrf' 



vN' 



C 



t + r] I \\ws\\ids 


t — 2 +7] \\wr\\idr \ / ||tfs||;^(is 



By Lemma 4.10, for any 77 > and p > 1 there exist positive constants C and 
7 so that for all A'' sufficiently large 

(4.13) E||Cf ir < C(l + llCo f )e''ll^"ll'e^^* = 2Ce''ll^"ll'e-''* . 

It remains to get control over the size of the perturbation v. Since v is adapted 
to the Wiener path. 



E 



vis)dWis] 



< / ^\v{s)fds<C / ^\Fsfds. 
Jo Jo 



Now since \\7riB{u,w)\\ < C\\u\\\\w\\ (see [EMSOl, Lemma A.4]), ||C/|| < 1 and 
llCf II = for t > 2, we see from (4.12) that there exists a C = C{N) such that 

for all s > 

1/2 



E||F,r <C l{,<2} + E||ti;,fE||C. 



By using (4.13) with p = 4, Lemma A.l from the appendix, to control E||ti;s||*^, 
and picking N sufficiently large, we obtain that for any rj > there is a constant 
C such that 

POD 

(4.14) E 



/•oo 

/ v{s)dWis) <Cexp{T]\\wof) . 
Jo 

Since (4.13) and (4.14) plug into (4.9), the result follows. 



D 
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4.6. The truly hypoelliptic setting: Proof of Proposition 4.3. We now 
turn to the truly hypoelliptic setting. Unlike the previous section, we allow 
for unstable directions which are not directly forced by the noise. However, 
Proposition 4.4 shows that the randomness can reach all of the unstable modes 
of interest, i.e. those in Ti. In order to show (4.1), we fix from now on ^ G W 
with ll^ll = 1 and we obtain bounds on (Wnfivj),^) that are independent 
of^. 

The basic structure of the argument is the same as in the preceding section 
on the essentially elliptic setting. We will construct an infinitesimal perturba- 
tion of the Wiener path over the time interval [0, t] to approximately match the 
effect on the solution wt of an infinitesimal perturbation of the initial condition 
in an arbitrary direction ^ & 7i. 

However, since not all of the unstable directions are in the range of Q, 
we can no longer infinitesimally correct the effect of the perturbation in the 
low mode space as we did in (4.12). We rather proceed in a way similar to 
the start of Section 4.3. However, since the Malliavin matrix is not invert- 
ible, we will regularize it and thus construct a v which compensates for the 
perturbation ^ only asymptotically as t -^ cx). Our construction produces a 
V which is not adapted to the Brownian filtration, which complicates a little 
bit the calculations analogous to (4.14). A more fundamental difficulty is that 
the Malliavin matrix is not invertible on any space which is easily identifiable 
or manageable, certainly not on Lq. Hence, the way of constructing v is not 
immediately obvious. 

The main idea for the construction of v is to work with a regularized 
version Mg^t = Mg^t + /3 of the Malliavin matrix Mg^t, for some very small 
parameter /3 to be determined later. The resulting M~^ will be an inverse "up 
to a scale" depending on (3. By this we mean that M~^ should not simply be 
thought of as an approximation of M^^ . It is an approximation with a very 
particular form. Theorem 4.12 which is taken from [MP06] shows that the 
eigenvectors with small eigenvalues are concentrated in the small scales with 
high probability. This means that M~^ is very close to M~^ on the large scales 
and very close to the identity times /3~^ on the small scales. Hence M~^ will 
be effective in controlling the large scales but, as we will see, something else 
will have to be done for the small scales. 

To be more precise, define for integer values of n the following objects: 

Jn = Jn,n+- i ^n = Jn+-,n+l ' -^n = An,n+- ' ^n ~ ^n^n i ^n = P + M„ . 

We will then work with a perturbation v which is given by on all intervals 
of the type [n + 2)'^ + 1]) ^-'^d by Vn G L?([n, n + 2],!^™') on the remaining 
intervals. 

We define the infinitesimal variation t;„ by 

(4.15) Vn = A*^M^^JnPn , 
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where we denote as before by pn the residual of the infinitesimal displacement 
at time n, due to the perturbation in the initial condition, which has not yet 
been compensated by v; i.e. pn = Jo,nC ~ ^o,n^o,n- From now on, with a slight 
abuse of notation we will write f „ for the perturbation of the Wiener path on 
[n, n + ^] and its extension (by 0) to the interval [n, n + I]. 

We claim that it follows from (4.15) that pn is given recursively by 

(4.16) pn+l = Jnf3M;;^JnPn , 

with po = C- To see the claim observe that (4.16) implies Jn,n+iPn = JnJnPn = 
JnAnVn + Pn+i- Using this and the definitions of the operators involved, we 
see that indeed 

N-l N~l 

Aq,NVo^N = 2_^ J {n+l),N JnAnVn = ^ {Jn,NPn — J(n+l),NPn+l) 
n=0 n=0 

= Jo,N^ — Pn ■ 

Thus, at time N, the infinitesimal variation in the Wiener path wq^at corre- 
sponds to the infinitesimal perturbation in the initial condition ^ up to an 
error p^v- 

It therefore remains to show that this choice of v has desirable properties. 
In particular we need to demonstrate properties similar to (4.13) and (4.14). 
The analogous statements are given by the next two propositions whose proofs 
will be the content of Sections 4.7 and 4.8. Both of these propositions rely 
heavily on the following theorem obtained in [MP06, Th. 6.2]. 

Theorem 4.12. Denote by M the Malliavin matrix over the time interval 
[0, 2] CLnd define Ti as above. For every a,r],p and every orthogonal projection 
iT£ on a finite number of Fourier modes, there exists C such that 

(4.17) P{{Mip,ip) < e\\ip\\l) < CeP eicp{7]\\wof) , 

holds for every {random) vector (p G 7i satisfying ||vr£(/9|| > aUv'lli almost surely, 
for every e G (0, 1), and for every wq € 7i. 

The next proposition shows that we can construct a v which has the 
desired effect of driving the error pt to zero as t ^ 00. 

Proposition 4.13. For any t] > 0, there exist constants /3 > and C > 
such that 

(4.18) E||p^||i°<^^"P^^"^°"') 



2N 

holds for every N > 0. [Note that by increasing (3 further, the 2^ in the 
denominator could be replaced by K for an arbitrary K > 2 without altering 
the value of C .) 
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However for the above result to be useful, the "cost" of shifting the noise 
by V (i.e. the norm of v in the Cameron-Martin space) must be finite. Since the 
time horizon is infinite, this is not a trivial requirement. In the "essentially 
elliptic" setting, it was demonstrated in (4.14). In the "truly hypoelliptic" 
setting, we obtain 

Proposition 4.14. For any t] > 0, there exists a constant C so that 

(4.19) E/ v^,,dW{s) < e''ll-«ll'V(E||/,„||iO)i. 

•^0 ^ n=0 

{Note that the power 10 in this expression is arbitrary and can he brought as 
close to 2 as one wishes.) 

Plugging these estimates into (4.9), we obtain Proposition 4.3. Note that 
even though Proposition 4.3 is sufficient for the present article, small modifi- 
cations of (4.9) produce the following stronger bound. 

Proposition 4.15. For every rj > and every 7 > 0, there exist con- 
stants Cfj^-y such that for every Frechet differentiable function ip from 7i to R, 

||VP„(^(u;)||<exp(r/||u;f)(c,,^^(P„|^|2)(t,)+7«^(p„||Vv5||2)(u;: 

for every w £ 7i and n G N. 

Proof. Applying Cauchy-Schwarz to the terms on the right-hand side of 
the penultimate line of (4.9) one obtains 

/ /■" 2 \ 1/2 / \ 1/2 

|(VP„99,e)|< (E / Vo,sdW{s) VrM^) + (BWpnW'^VnWV^f) . 

It now suffices to use the bounds from the above propositions and to note that 
the right-hand side is independent of the choice of ^ provided ||^|| = 1. D 

4.7. Controlling the error: Proof of Proposition 4.13. Before proving 
Proposition 4.13, we state the following lemma, which summarizes the effect 
of our control on the perturbation and will be proved at the end of this section. 

Lemma 4.16. For every two constants 7,7/ > and every p > 1, there 
exists a constant Pq > such that 



E(||/>„+i|n.F„)<7e''""'"" \\p: 



'n\ 



IP 



holds almost surely whenever P < Pq. 
Proof of Proposition 4.13. Define 



lU l|10 
^ _ llPn+lll 

" ~ II inn 

1 1 r n 1 1 



ERGODICITY OF THE 2D NAVIER-STOKES EQUATIONS 1021 

with the convention that C„ = if /?„ = 0. Note that since \\po\\ = 1, one has 
II/Oa^ ||^° = nn=o ^n- We begin by establishing some properties of C„ and then 
use them to prove the proposition. 

Note that ||/3M~^|| < 1 and so, by (4.8) and (4.16), for every rj > there 
exists a constant C^ > such that 
(4.20) 

C < W 1 RA/r~'^ T W'^^ < II 1 ll^f^ll 7" ll^*^ < PM-nl n / Hid II 2 We 4- r* 

^n ^ W'JnP-^^J-n '^nW ^|Kn|| ||<^n|| ^ cXp I f/ / || u;^ 111 "'' i^ ^-^r; 

Jn 

almost surely. Note that this bound is independent of /3. Next, for given values 
of Tj and R > 0, we define 

r e-'?^ if \\wnf > 2R, 

Cn,R = \ „ 

[ e^ Cn otherwise. 

Obviously both C„ and Cn,R are .7-"„+i-measurable. Lemma 4.16 shows that 
for every R > r]~^, one can find a /? > such that 

(4.21) ^{Cl^R I -^n) < ^ , almost surely . 

Note now that (4.20) and the definition of Cn^ji immediately imply that 

(4.22) Cn<Cn,Rexp\r] \\ws\\ids + r]\\wn\\'^ + Cr, - r]R] , 

almost surely. This in turn implies that 

Af-l A-1 7V-1 ^n+1 



llCn<ll CIr + n e^P(2^ / W^sWlds + 21^\\Wnf + 2Cr, " 2r/i? 

N-1 N-1 

n=0 ra=0 

+ expU7? / ll^slli (is + 2iV((7^ -7?i?) 



Now fix ry > (not too large). In light of (4.6) and (4.7), we can then choose R 
sufficiently large so that the last two terms satisfy the required bounds. Then, 
we choose /3 sufficiently small so that (4.21) holds and the estimate follows. D 

To prove Lemma 4.16, we will use the following two lemmas. The first is 
simply a consequence of the dissipative nature of the equation. Because of the 
Laplacian, the small scale perturbations are strongly damped. 

Lemma 4.17. For every p > 1, every T > 0, and every two constants 
7,r/ > 0, there exists an orthogonal projector tt^ onto a finite number of Fourier 
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modes such that 

(4.23) E||(l - 7r,)Jo,rr < jexp{7]\\wof) , 

(4.24) E|| Jo,r(l - vr,)r < 7exp(??||w;of ) • 

The proof of the above lemma is postponed to the appendix. The second 
lemma is central to the hypoelliptic results in this paper. It is the analog of 
(4.14) from the essentially elliptic setting and provides the key to controlling 
the "low modes" when they are not directly forced and Girsanov's theorem 
cannot be used directly. This result makes use of the results in [MP06] which 
contains the heart of the analysis of the structure of the Malliavin matrix for 
equation (2.1) in the hypoelliptic setting. 

Lemma 4.18. Fix ^ G TY and define 

Then, for every two constants 7, r/ > and every low-mode orthogonal projec- 
tor TTi, there exists a constant /3 > such that 

Remark 4.19. Since one has obviously that ||C|| < ||<^o?||; this lemma tells 
us that applying the operator /?(/? + Mq)^^ (with a very small value of (3) to a 
vector in Ti. either reduces its norm drastically or transfers most of its "mass" 
into the high modes (where the cutoff between "high" and "low" modes is 
arbitrary but influences the possible choices of /?). This explains why the 
control V is set to for half of the time in Section 4.6. In order to ensure that 
the norm of pn gets really reduced after one step, we choose the control in 
such a way that /3(/3 + Mn)~^Jn is composed by J„, using the fact embodied 
in Lemma 4.17 that the Jacobian will contract the high modes before the low 
modes start to grow out of control. 

Proof of Lemma 4.18. For a > 0, let A^ denote the event ||vr^Cll > Q^HCIIi- 
We also define the random vectors 

where lj is the chance variable and XA is the characteristic function of a set A. 
It is clear that 

E||7r,Cair<aPE||C||?. 

Using the bounds (4.8) and (4.6) on the Jacobian and the fact that Mq is a 
bounded operator from T^i (the Sobolev space of functions with square inte- 
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grable derivatives) into 7ii, we get 

(4.25) EllvrKair < «^E||C||? < a^E||Joe||? < ^e^^^^^^^'UV , 

(with Tj and 7 as in the statement of the proposition) for sufficiently small a. 
From now on, we fix a such that (4.25) holds. One has the chain of inequalities 

(4.26) (Ca, MoCa) < (C, MoC) < (C, (Mo + m 



From Theorem 4.12, we see furthermore that, for every po > 0) there exists a 
constant C such that 

P{{MoCa,Ca) < eWCaWl) < CsP" e^p{v\\wof) , 

holds for every wq £ 7i and every e £ (0, 1). Consequently, 

/ \\/- ||2 1 s 

P(p^ >-)< P((MoCa,Ca) < ePWUJ) < C^P"eP^'eMv\\wof) , 

where we made use of (4.26) to get the first inequality. This implies that, for 
every p,q > 1, there exists a constant C such that 

(4.27) Ef^^) < C^''exp(i]\\wof) . 



Since Hvt^Cq-H < ||Ca||i and 

E||ail?<WE 



V hMW'p^ 



it follows from (4.27) and the bound (4.8) on the Jacobian that, by choosing 
P sufficiently small, one gets 

(4.28) EiKKair < ^e^ii-^'Hieir • 

Note that E||7r£("||^ = E||7r^(^air + E||7r^(^Q,||*' since only one of the previous two 
terms is nonzero for any given realization uj. The claim thus follows from (4.25) 
and (4.28). D 

Using Lemma 4.17 and Lemma 4.18, we now give the 

Proof of Lemma 4.16. Define Cn = (3M~^JnPn, so that Pn+i = JnCn- 
It follows from the definition of M„ and the bounds (4.8) and (4.6) on the 
Jacobian that there exists a constant C such that 



E(||C„in-^n)<Ceil 



'""ll lln IIP 
1 1 Pn 1 1 
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uniformly in /3 > 0. Applying (4.24) to this bound yields the existence of a 
projector tt^ on a finite number of Fourier modes such that 



E(|| J„(l - 7r,)C„in -^n) < Te^'ll'-il'llp. 



IP 



Furthermore, Lemma 4.18 shows that, for an arbitrarily small value 7, one can 
choose f3 sufficiently small so that 

E(||7r,Cn|ri-^n)<7e^ll"'"iniPnir- 

Applying again the a priori estimates (4.8) and (4.6) on the Jacobian, we see 
that one can choose 7 (and thus /3) sufficiently small so that 

E(||J„7r,C„in-^n)<7e''"""'niPn|r, 
and the result follows. D 

4.8. Cost of the control: Proof of Proposition AAA. Since the process fo,s 
is not adapted to the Wiener process VF(s), the integral must be taken to be a 
Skorokhod integral. We denote by T)sF the Malliavin derivative of a random 
variable F at time s (see [Nua95] for definitions). Suppressing the dependence 
on the initial condition tt;, we obtain from the definition of the Skorokhod 
integral and from the corresponding Ito isometry (see e.g. [Nua95, p. 39]) 

-No ^ /-n+i /.n+i 



E 



/ v{s)dW{s) <E||vo,7vf + J^ / / mVsVn{t)fdsdt. 

JO „_Q^ra Jn 



(Remember that Vn{t) = on [n + 2 , Ji- + 1]-) In this expression, the norm ||| • ||| 
denotes the Hilbert-Schmidt norm on m x m matrices, so that one has 

/ / B\l'D,Vn{t)fdsdt = ^ EWVlvJ^ds, 

Jn Jn -I Jn 

where the norm ||-|| is in l?{[n,n + ^],'R'"'') and !D* denotes the Malliavin 
derivative with respect to the i component of the noise at time s. 

In order to obtain an explicit expression for T)\v^, we start by computing 
separately the Malliavin derivatives of the various expressions that enter into 
its construction. Recall from [Nua95] that T)\wt = Js,tQ^i for s < t. It follows 
from this and the expression (4.2) for the Jacobian that the Malliavin derivative 
of Js,tC is given by 

dt1)iJs,tC = i^A1)Us4 + B{wt, 1)iJs,tO + B{Jr,tQei, Js4) . 

From the variation of constants formula and the expression (4.4) for the process 
K, we get 

f Kr,tiQei,Js,rO ifr>s, 
(4.29) 'DrJs^ = < 

I Ks,t{Jr,sQei-,i) if r < s. 
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In the remainder of this section, we will use the convention that if A : TCi -^ 7i2 
is a random linear map between two Hilbert spaces, we denote by T)\A : 7ii -^ 
0.2 the random linear map defined by 

{VlA)h = {Vs{Ah),e^) . 

With this convention, (4.29) immediately yields 

(4.30) TilJj^w = Krn+^{Jn,rW,Qei) for r G [n,n + -]. 

Similarly, we see from (4.29) and the definition of ^4^ that the map "D^A^ given 

by 

j'r 

(4.31) ■I)iA^h= K^^^^i_{Js,rQh{s),Qei))ds 

J n 

+ / -f^r,n+i {Qh{s), Jr,sQei)) ds . 

We denote its adjoint by D^yl* . Since M„ = (3 + AnA^, we get from the chain 
rule 

•DIM-' = -M-'({^lA,,)Al + A^{VIAI))m-' . 



Since pn is .7>j-measurable, one has !D^p„ = for r > n. Therefore, combining 
the above expressions with the Leibniz rule applied to the definition (4.15) of 
Vn yields 

D>„ = (du:)m„"' JnPn + ^:m-^(!d: j„)p„ 

-A*^M-'({VlA^)Al + A^{'DlA*^))M-UnPn. 



Since M„ = /3 + A^A'^, one has the almost sure bounds 

II 4* M^-*-/^!! < 1 WM"^/^ A II < 1 llM^-^/^ll < R~^/^ 

This immediately yields 

111)*?; II < 3/9^^117)* A Ml WWn II -I- /l^^/^IIT)* 1 Mr) II 
ll-i-'s^nll — "J/^ ll-^s^nllll'^nllllp'nll T^ P II -^s'^nll llPn II • 

Combining this with (4.31), (4.30), and Lemma 4.10, we obtain, for every 
r/ > 0, the existence of a constant C such that 



E||D>„f < Ce'?""'" /3-^(E||p„ 



Applying Lemma 4.10 to the definition of Vn we easily get a similar bound for 
E||w„|p, which then implies the quoted result. 
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5. Discussion and conclusion 

Even though the results obtained in this work are relatively complete, 
they still leave a few questions open. 

Do the transition probabilities for (2.1) converge towards the invariant 
measure and at what rate? In other words, do the solutions to (2.1) have 
the mixing property? We expect this to be the case and plan to answer this 
question in a subsequent publication. 

What happens iffi^TC and one starts the system with an initial condition 
Wo £ 7i\ Ti.7 If the viscosity is sufficiently large, we know that the component 
of wt orthogonal to 7i will decrease exponentially with time. This is however 
not expected to be the case when i/ is small. In this case, we expect to have (at 
least) one invariant measure associated to every (closed) subspace V invariant 
under the flow. 

Appendix A. A priori estimates for the Navier-Stokes equations 

Note. The letter C denotes generic constants whose value can change 
from one line to the next even within the same equation. The possible depen- 
dence of C on the parameters of (2.1) should be clear from the context. 

We define for a G R and for w a smooth function on [0, 27r]^ with mean 
the norm \\w\\a by 

where of course Wk denotes the Fourier mode with wavenumber k. Define fur- 
thermore {K,w)k = —iwkk^/WkW^, B{u,v) = {u-'V)v and 5 = {s = (si, S2, S3) G 
KI : J2 Si > 1, s j^ (1, 0, 0), (0, 1, 0), (0, 0, 1)}. Then the following relations are 
useful (cf. [CF88]): 

(A.l) {B{u,v),w) = -{B{u,w),v) ifV-u = 0, 

(A.2) \{B{u,v),w)\ <C\\u\\s,\\v\\i+s2\\w\\s3, (81,82,83) £ S, 

(A. 3) ||A^^t||a = ||'u||q-1, 

(A.4) ||u||| < e||ti||^ + e"^^||ii||^ if < a < /3 < 7 and e > 0. 

Before we turn to the proof of Lemma 4.10, we give the following essential 
bound on the solutions of (2.1). 

Lemma A.l. There exist constants rjQ > and C > 0, such that for every 
t > and every rj £ (0, ryo], the bound 

(A.5) Bexp{r]\\wtf) < Cexp(r/e"''*||ti;of ) 

holds. 
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Proof. From (A.l) and Itos formula, we obtain 

\\wt\f - \\wo\\'^ + 2u \\wr\\ldr= {wr , Q dW (r)) + Sot , 
Jo Jo 

where Sq = tr QQ*. Using the fact that Iju'rlli > Uti^rpj we get 

^ Jo Jo 

There exists a constant a > such that z^Htf^lp > ^IIQ*^^!^, so that [Mat02a, 
Lemma A.l] imphes 

V u a / 

Note now that if a random variable X satisfies P{X > C) < 1/C^ for all 
C > 0, then EX < 2. The bound (A. 5) thus follows immediately with for 
example tjq = a/2 and C = 2exp(^|i). D 

We now turn to the proof of Lemma 4.10. 

Proof of Lemma 4.10. Point 1. From (A.l) and Ito's formula, for any 
r/ > we obtain 

^11^*11 + V^ \\wr\\i dr — ri£o{t — s) 

= riW-WsW"^ + rj {wr , Q dW {r)) — rju \\wr\\idr 

J s J s 

where Sq = tr QQ*. Denote by M{s,t) the first two terms on the right-hand 
side of the last expression and set N{.s,t) = M{s,t) — rju J^ WwrWldr. Now 
observe that with a as in the proof of Lemma A.l above, one has N{s,t) < 
M(s, t) — ^{M)(s, t) where {M)(s, t) is the quadratic variation of the continu- 
ous L^-martingale M. Hence by the standard exponential martingale estimate, 
P(supi>^iV(s,t) > kite's) < e^p{r]\\wsf-^) for ah s > 0. Here we use the 
notation J^g to denote the filtration generated by the noise up to the time s. 
Thus, for all rj £ (0, a/2] and s > 0, 

E expiry sup ( ||tt;( II + '^ / \\wr\\idr — £(){t — s)] J^g) < 2exp(77||ti;s|| )• 

Choosing r/o as above and using Lemma A.l to bound the expected value of 
the right-hand side, we complete the proof. 

Point 2. Taking conditional expectations with respect to J-'n-i on the 
left hand side of (4.7) and applying Lemma A.l, one has 

N N-1 

Bexp(ri^\\wn\\'^j < CBexpUie'" \\wN -iW"^ + V^ \\wnf 

n=0 n=0 
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Applying this procedure repeatedly, one obtains 

N 



Eexp(r?^||it)„|H < C^exp(ary||ti;o||^) , 



n=0 

where a = ^^=q e~^"'. This computation is vahd, provided ar] is smaher than 
r/o, so the result fohows by taking rji = rjo/a. 

Point 3. We define ^t = ■Jo,tS,o for some .^o £ W. The evolution of ^j is 
then given by (4.2). We thus have for the H-norm of S, the equation 

dtUtf = -2iy\\V^tf + 2{B{)C^t,wt),^t) ■ 

Equation (A.l) yields the existence of a constant C such that 2| {B{ICh, w),C) \ ^ 
C||w||i||^||||Clli/2 foi' example. By interpolation, we get 

(A.6) 2| {B{ICh,w),0 I < HICII? + §-\\Cf + hwWlWhf , 

and therefore 

(A.7) dMtf < -^l|V6f + -^Utf + hwtWlUtf , 

for every ry > 0. This yields (4.8). 

Point 4. This bound follows in a rather straightforward way from (A. 8) 
which is in the next proof. Standard Sobolev estimates and interpolation 
inequalities give for the symmetrized nonlinearity B the bound 

\\B{u,w)\\ < C{\\u\\ii2\\w\\i + ||M||i||tt;||i/2) 

<C(||ur/2||n||}/'||u;||i + ||u;r/2||u;||}/2||u||i). 

Combining this with the definition (4.4) of Kg^ and bound (A. 8) yields for 

s,tG[0,l] 



\Ks,t\\<C f \\Jr,t\\\\JsAT\\JsA\'^^dr 
<Cexp(r/ / ||ti;r||idr 



where we used the integrability of |r — s| ^' "^ in the second step. This concludes 
the proof of Lemma 4.10. D 

Proof of Lemma 4.17. In order to get (4.23), we show that with the above 
notation, one can get bounds on ||^t||i as well. To achieve this we define, for 
a constant e > to be fixed later, Ci = llCilP + ^^ll'^tlli- Using (A.7) to bound 
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the derivative of the first term and combining (A. 2) with (A. 3) for the other 
terms, we then get in a straightforward way 



dtCt <ie- ^)||V6f + ^||6f + IWwtWlUtl 

^t||l/2 • 



-2teiy\Ml + 2teC\\wt\\iUthUt 



By (A. 4), we get 



2cihiiiiicii2iiciii/2<2v^ikiiiiieii2ii^iii + |;ikiiiiieii2iieii 

This immediately yields 

dtCt <{s- u)\Ml + -^Utf + (? + 4^) WMlUtf 

-teu\\it\\l + ter^\\wt\\l\M\. 
If we take e sufficiently small (of the order "rfv^), we get 

dtQt<[-^+Tl\\wt\\t)Qt. forte [0,1], 
and therefore 

(A.8) \\Mo\\l<jew{vJ \\ws\\lds)Uof, 

for some (possibly rather large) constant C. If we now define ttn as the or- 
thogonal projection on the set of Fourier modes with \k\ > N, 

hN^tW < ^||6lli • 

The bound (4.23) immediately follows by taking vr^ = 1 — tt^v for N sufficiently 
large. 

We now turn to the proof of the bound (4.24). We define ire as above (but 
reserve the right to choose the precise value of N later) and set Q = iTi^t and 
^^ = (1 — 7re)^t- With this notation, (4.24) amounts to obtaining bounds on 
||.^t|| with ^Q = 0. Using the identity (A.l) we have 

dM^f = -2,yUX-2{BilC^t,^^),wt)-2{BilCwt,^^),^t. 
Applying (A. 2) to the right-hand side allows us to get the bound 

Stiidf < -2Hicf II? + 2 (i?(/ce,^ ^i), ef ) + cii«;tiii/2iief iiiiie'ii , 
9tiie*'f<-2HiCt'ii? + cik*iii/2iie.'iiiii6ii. 
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We then bound the first fine using (A. 6) and the second fine using ||^^||i > 
N'^UtW^- We thus obtain 

(A.9) mff < (^ + vwwtwi) uff + c\\w,\\i^,un' , 

for an arbitrary value of r] and for a constant C depending on u but independent 
of N and rj. Using the a priori bound from point 3 above for the Jacobian ^ 
and the interpolation inequality ||tL's||^/2 < ||tt)s||||ti's||i immediately produces 
the bound 

t 



Jo 



Hu.hu2 , C{T)\\^f^\\ ft 11^^ ||.^^ 



<e-^^1egr+ T"" e 



N \i ./n ^(=\n t] 



/ ll^sllf ds sup \\Ws 
Jo se\o,t] 



<||^o^||2(e-^^* + ^exp(ry' f \\ws\\Us) sup \\w. 
^ ^^ ^ Jo ' sefo.tl 



' se[o,t] 

for an arbitrary r/' > r]. Combining this with the bound of point 1 above shows 
that, for every ry, every 7, every p, and every T, there exists a constant Nq 
such that 

for all t G [0,T] and all N > Nq. Since Co ~ ^ ^J assumption, it follows from 
(A.9) that 

\\df<C f\w{^^ + v[\\wr\\ldr)\\ws\\l/Mfds 

"/ U / ^ s 



' / ||ttJs||^ ds sup llt^sll . 
Jo se\o,t] 



m 

The required bound (4.24) now follows easily by taking expectations and using 
the previous bounds. D 
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